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ABSTRACT 

We report the detection of luminous C0(J=2^-1), C0(J=3— ^2), and C0(J=4— s>3) emission in the 
strongly lensed high-redshift quasars B1938+666 (z=2.059), HE 0230-2130 (z=2.166), HE 1104-1805 
(z=2.322), and B1359+154 (z=3.240), using the Combined Array for Research in Millimeter-wave 
Astronomy. 61938+666 was identified in a 'blind' CO redshift search, demonstrating the feasi- 
bility of such investigations with millimeter interferometers. These galaxies are lensing-amplified 
by factors of /xl— 11-170, and thus allow us to probe molecular gas in intrinsically fainter galaxies 
than currently possible without the aid of gravitational lensing. We report lensing-corrected intrinsic 
CO line luminosities of L^q=0. 65-21 x 10^ K kms^^pc^, translating to H2 masses of M(H2)=0.52- 
17x10^ (aco/O-8)M0. To investigate whether or not the AGN in luminous quasars substantially 
contribute to Lfir, we study the L^q-Lfir, relation for quasars relative to galaxies without a lumi- 
nous AGN as a function of redshift. We find no substantial differences between submillimeter galaxies 
and high-z quasars, but marginal evidence for an excess in LpiR in nearby low-LpiR AGN galaxies. 
This may suggest that an AGN contribution to LpiR is significant in systems with relatively low 
gas and dust content, but only minor in the most far- infrared- luminous galaxies (in which LpiR is 
dominated by star formation). 

Subject headings: galaxies: active — galaxies: starburst — galaxies: formation — galaxies: high- 
redshift — cosmology: observations — radio lines: galaxies 



1. INTRODUCTION 

Molecular line emission from galaxies at high redshift 
(z>2) has proven to be a unique tool to investigate the 
stellar mass buildup in the often heavily dust-enshrouded 
star-forming environments of massive galaxies at early 
cosmic times. Molecular gas (CO) was detected in 26 
2>2 quasar host galaxies to dat e (see, e.g., review by 
Solomon & Vanden Bout l2005| ). all of which are lu- 
minous in the rest-frame far infrared (FIR; typically 
iFiR>10"Lo; e.g., Wang et al. 120081 ). These stud- 
ies typically reveal molecular gas masses of few times 
lO^°A'/0, indicating gas-rich host galaxies for the lu- 
minous active galactic nuclei (AGN) that are observed 
at optical wavelengths. The high FIR luminosities are 
thought to be dominated by dust heating from young 
stars formed in intense starburst events comparable to 
those found in submillimeter galaxies (SMGs), with star 
formation rates (SFRs) exceeding lOOOMQyr"^. These 
starburst are fueled by the massive molecular gas reser- 
voirs, wh ich ca n maintain them for >10^yr (e.g., Riech- 
ers et al. 120081 ). The connection between gas mass and 
star formation is reflected in the relation betw een CO 
and FIR lumin osities (e.g., Sanders et al. llQQlt Gao & 
Solomon [200l . which in turn can be used to constrain 
the contributio n of t he AGN to Lpir in quasars (e.g., 
Riechers et al. 120061 ). The presence of both luminous 
AGN and massive starbursts make these objects the ideal 
candidates to better understand the AGN-starburst con- 
nection back to early cosmic times. 
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The currently most efficient way to investigate the 
molecular gas properties of quasar host galaxies with 
somewhat less extreme star formation events is by facil- 
itating the flux magniflcation provided by gravitational 
lensing. Besides the flux boost, the spatial magnification 
of lensed quasars also aids in resolving these systems, al- 
lowing us to probe the gas reservoirs at smaller physical 
scales. About one third of the CO-detected z>2 quasar 
host galaxies are gravitationally lensed (9 sources) , with 
typical magnification fac tors b etween a few and '-~^30 
(Solomon & Vanden Bout 120051). To study the gas prop- 
erties of high redshift AGN host galaxies in more detail, 
we here enhance this sample by four galaxies with magni- 
fication factors of 10.8-173, including double, quadruple, 
sextuple, and Einstein ring lens configurations. 

In this paper, we report the detection of C0(J=2— >1), 
CO(J=3-^2), and CO(J=4^3) emission in the far- 
infrared-luminous, strongly lensed quasar host galax- 
ies of B1938+666 (z=2.059), HE0230-2130 (z=2.166), 
HE 1104-1805 (z=2.322), and B1359+154 (z=3.240), us- 
ing the Combined Array for Research in Millimeter- 
wave Astronomy (CARMA). We use a concor- 
dance, flat ACDM cosmology throughout, with 
gn= 71 kms~^M pc~\ r2M=0.27, and rJA=0.73 (Spergel 
etal. [2003l[2007l) . 

2. OBSERVATIONS 

2.1. Targeted CO Searches 

We used CARMA to observe the CO(J=3^2) 
transition line (z/rost =345. 7960 GHz) toward the z>2 
quasars HE 0230-2130, HE 1104-1805, and B1359-hl54, 
and the CO(J=4-^3) line (1/^081=461.0408 GHz) toward 
B1359-I-154, all of which are redshifted to the 3 mm at- 
mospheric window (see Tab.[l]for redshifted frequencies). 
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TABLE 1 
Summary of observations. 



target 


line 


^opt 


[GHz] 


sideband 


configuration 


dates 


tracks 


calibrator 


^on/^tot 

[hr] 


B1938+666 


CO(J=2-s>l) 




75.3638 




C/D/E 


2010 Jul 05-Oct 10 


3/10/4 


1849+670 


43.6/69.9 




CO(J=3-s>2) 




113.0425 














HE 0230-2130 


CO(J=3-s>2) 


2.163 


109.360 


USB 


D 


2009 Aug 01-Sep 02 


5 


0132-169 


13.3/22.9 


















0204-170 




HE 1104-1805 


CO(J=3-s>2) 


2.319 


104.155 


USB 


C 


2010 Mar 21/22 


2 


1127-189 


6.6/10.6 


B1359+154 


CO(J=3-s>2) 


3.235 


81.652 


LSB 


D/E 


2009 Jul 07-Aug 15 


1/2 


1357+193 


5.9/10.6 




CO(J=4-s>3) 




108.864 


USB 


D/E 


2009 Jun 29-Jul 27 


2/2 




7.2/13.3 


total 














31 




76.6/127.3 



Tuning frequency (where apphcablc), corresponding to zero velocity in Figs. [Tl andFf] 
Time on source/total. 



All targets were observed with 14 or 15 antennas (cor- 
responding to 91 or 105 baselines per antenna configu- 
ration) for a total of 14 tracks in the C, D, and E con- 
figurations between 2009 June 29 and 2010 March 22, 
amounting to a total observing time of 57 hr (33 hr on 
source). 

Weather conditions scaled between acceptable and 
excellent for observations at 3 mm wavelengths. The 
nearby quasars J0132-169, J0204-170, J1127-189, and 
J1357-I-193 were observed every 15 minutes for secondary 
amplitude and phase calibration. The strong calibra- 
tor sources J0423-013, J1058-h015, 3C446, and 3C273 
were observed at least once per track for bandpass and 
secondary flux calibration. Absolute fluxes were boot- 
strapped relative to Mars, Uranus, MWC349, or 3C84 
(when no planet was available) . Pointing was performed 
at least every 2-4 hr on nearby sources, using both ra- 
dio and optical modes. The resulting total calibration is 
estimated to be accurate within ^^15%. 

All observations were carried out with the previ- 
ous generation correlator. The CO(J=3-^2) line in 
B1359-I-154 was centered in the lower sideband (LSB). 
All other lines were centered in the upper sideband 
(USB), at intermediate frequencies of 2.5 GHz. Three 
bands with 15 channels of 31.25 MHz (-86-115 kms'^) 
width each were centered on the tuning frequencies. 
The bands were overlapped by 2 channels to improve 
calibration of the correlated dataset, leading to an ef- 
fective bandwidth of 1281.25MHz (-3500-4700 kms'^) 
per sideband. 

For data reduction and analysis, the MIRIAD package 
was used. The final plots were created with the GILDAS 
package. All data were imaged using 'natural' weighting, 
yielding synthesized beam sizes of 7.3"x4.2", 2.5"xl.9", 
10.1"x7.8", and 9.4"x6.8" for the CO(J=3^2) observa- 
tions in HE 0230-2130, HE 1104-1805, and B1359-H154 
and the CO(J=4-)-3) observations of B1359+154, re- 
spectively. The final rms noise values are 0.52, 0.65, 
0.97, and 0.71 mJy beam"! over 1114, 450, 344, and 
430 kms-i (406.25, 156.25, 93.75, and 156.25MHz). Av- 
eraging over all line-free data (LSB-I-USB) in the obser- 
vations of B1359-I-154 (both line setups) yields a beam 
size of 9.8"x7.3" and an rms noise of 0.12 mJy beam~^. 

2.2. 'Blind' CO Search 

We also used CARMA to observe the very 
submillimeter-bright, radio-loud quasar B1938+666, for 
which the redshift was previously unknown. Thus, 
we searched the entire frequency range from 82.4 to 



115.3 GHz for CO line emission, facilitating the large 
bandwidth of 3708.096 MHz per sideband (eight bands 
with 95 channels of 5.208 MHz with per sideband, over- 
lapping bands by typically 6 channels to reduce sensitiv- 
ity losses due to bandpass rolloff) of the new correlator. 
Observations were carried out for 14 tracks in the D and 
E configurations between 2010 July 05 and September 
08, amounting to a total observing time of 58 hr (37 hr 
on source). A total of 6 (partially overlapping) frequency 
settings were used to cover the 82.4-115.3 GHz frequency 
range. Observations were typically carried out in pairs 
of frequency settings with an IF frequency of 3.6 GHz, 
where the LSB of the second setup was used to fill the 
gap between sidebands of the first setup (and, vice versa, 
the USB of the first setup fills the IF gap of the second 
setup) to achieve continuous frequency coverage. 

After identification of a CO line, a seventh frequency 
setup was used to observe a second CO line, redshifted 
to ~75.4GHz, for confirmation. This frequency is almost 
10 GHz below the nominal tuning range of CARMA's 
3mm receivers (85-116 GHz). We thus placed the line 
in the LSB, using the extended IF range with an IF fre- 
quency of 7.5 GHz (same bandwidth as above), and plac- 
ing the local oscillator (LO) at 82.874 GHz0 With this 
setup, we successfully extended the covered frequency 
range down to —74.1 GHz, with decreasing sensitivity 
below —76 GHz. These observations were carried out for 
3 tracks in C configuration between 2010 September 27 
and October 10, amounting to a total observing time of 
12 hr (6.6 hr on source). 

Weather conditions scaled between acceptable and ex- 
cellent for observations at 3 mm wavelengths. The nearby 
quasar J1849+670 was observed every 12-20 minutes 
for secondary amplitude and phase calibration. The 
strong calibrator sources J1751-h096, J2015-t-372, 3C273, 
3C345 and 3C454.3 were observed at least once per track 
for bandpass and secondary fiux calibration. Absolute 
fluxes were bootstrapped relative to Neptune and Mars. 
Pointing was performed at least every 2-4 hr on nearby 
sources, using both radio and optical modes. The result- 
ing total calibration is estimated to be accurate within 
-15% (20% for setup 7). 

For data reduction and analysis, the MIRIAD pack- 
age was used. The final plots were created with 
the GILDAS package. All data were imaged using 
'natural' weighting, yielding synthesized beam sizes 

^ The USB falls into the frequency range covered by the previous 
setups. 
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TABLE 2 
Observed Line Parameters. 



target line zco Sv AijpwHM Ico ^J+l,j'' 

[mjy] [kms^-'-] [Jy kms^-'-] 

B1938+666 CO(J=2-s>l) 2.0590±0.0003 9.9±2.6 366±124 3.8±1.1 1.0±0.3 
CO(J=3^-2) 16.3±1.8 529±75 9.1±1.1 
HCN(J=3^.2) (3.1±1.1) (496±192)'' (1.62±0.54)* 
HCO+(J=3^2) (1.6±1.1) (496±192)'' (0.86±0.45)* 
HNC(J=3^2) (1.9±1.1) (496±192)'' (1.00±0.46)* 
C2H(Af=3-i.2) (1.7±1.1) (496±192)b (0.89±0.46)* 
CN{Ar=3^-2) (1.5±1.1) (1.0±0.7)* 
HE 0230-2130 CO(J=3^-2) 2.1664±0.0005 11.1±1.5 705±123 8.3±1.2 
HE 1104-1805 CO(J=3^2) 2.3221±0.0004 16.1±2.4 441±81 7.5±1.2 
B1359+154 CO(J=3^-2) 3.2399±0.0003 5.6±1.7 198±92 1.2±0.4 1.2±0.5 
CO(J=4-»3) 10.0±1.6 237ifc47 2.5ifc0.4 

Note. — *: Considered not detected. Corresponds to 3cr upper limits of 1.6 (HCN), 1.4 (HCO+, HNC, C2H), and 2.1 Jy kms~^ (CN). 
^ Line brightness temperature ratio CO( J— n+1— >n)/CO(t/— n— >n— 1), where n— 2/3 (r— 1: thermalized) . 
^ Fixed to a common linewidth. The width of the C2H line is corrected for broadening due to hypcrfine structure. 



of 2.8"x2.3" and 4.5"x3.9" at 75.4 and 113.0GHz, 
and 4.5"x4.0" over the full, 41.2 GHz wide band- 
pass (with a small gap around 80 GHz, and exclud- 
ing line emission). Due to varying effective expo- 
sure time, atmospheric conditions and receiver noise, 
the sensitivity varies somewhat over the full spectral 
bandpass, in particular toward the edges of the 3 mm 
band. The final rms noise values are 24^Jybeam~^ 
over the full bandpass, 0.77mJy beam~^ over 203 MHz 
(808 kms-i) at 75.4GHz, and 0.52 mJy beam"! over 
292MHz (773 kms"!) at 113.0GHz. 

3. RESULTS 
3.1. B 1938+666 
3.1.1. Previous Results 

B1938+666 is a radio-selected gravitational lens with 
a small image separation of only 0.95" in the radio con- 
tinuum emission associated with its radio - loud active 
galactic nucleus (AGN; Patnaik et al. ll992D . The bright 
radio continuum emission (0.577±0.017 Jy at 1.4 GHz) 
is lensed into at least 3 components that break up into 
7 subcomponents along a partial Einstein ri ng, so me of 
which show polarized emission (King et al. Il997h . The 
radio components show a relatively steep spectral in- 
dex that varies between the different subcomponents 
(af gi2=-(l-0-0.5) from 1.612 to 5 GHz), suggesting a 
superposition of flat and steep spectrum emission asso- 
ciated with lensed core and jet emission. Its 1.6 /im con- 
tinuum emission is lensed into a full, almost perfectly 
symmetric Einstein ring with 0.95" diameter, but only 
the foreground lensing galaxy is detected at i band and 
shorter wavelengths (King et al. 1998). This shows that 
B1938-I-666 has properties that are clearly different from 
optically-selected, gravitationally lensed quasars. The 
lensing galaxy has a spectroscopic redshift of z=0.881, 
suggesting z>l.l for B1938-I-666 based on the lensing 
configuratio n, wit h a most likely redshift of 2=2.8 (Tonry 
& Kochanek|2003). Until now, the redshift of B1938+666 
was highly uncertain. In addition to its radio-loud 
AGN, B1938-I-666 also hosts a luminous dust bump, 
with S'45o^m=126±22mJy, 5'850Mm=34.6±2.0mJy, and 
S*! 3^jjj=14.7±2.0mJy (Barvainis & Ivison il002). Cor- 
recting for synchrotron emission from the AGN, Barvai- 
nis & Ivison (,2002.) estimate that 32.0±6.7mJy of the 



850 /im emission are due to dust heating, likely domi- 
nated by star formation. They also find a model-based 
lensing magnification factor of /zl=173. Its extreme sub- 
millimeter brightness, dust obscuration in the optical and 
small Einstein radius (making optical spectroscopic red- 
shift searches unfeasible) make B1938-)-666 an ideal can- 
didate to attempt a 'blind' CO redshift search. 

3.1.2. New Results: 'Blind' CO Redshift and Line Emission 

Scanning the whole 3 mm band from 82.4-115.3 GHz, 
we have identified a single CO line in B1938-|~666 at 
113.0425 GHz. The lack of detecting a second, lower- 
J CO line at 82.4<z<113.0GHz rules out any candi- 
date redshifts at z>3.0. Assuming a lower limit of 
z=1.7 as suggested by the lensing configuration, this 
leaves C0(J=3— ?>2) emission at z=2.059 as the only 
likely explanation for the detected line emission. We con- 
firmed this conclusion through a subsequent detection of 
C0(J=2^1) emission at 75.3638 GHz. 

In Fig. [1] the full covered spectral range is shown, 
along with higher spectral resolution, zoomed-in spectra 
of the emission lines, and emission line maps (contin- 
uum not subtracted). As both redshifted CO lines are 
close to the edges of the 3 mm band, the noise across the 
zoomed-in spectra is not flat, and substantially increases 
below ^75 GHz and above ^114 GHz. From Gaussian 
fitting to the line profiles, we derive line peak flux 
densities of 5'^=9.9±2.6 and 16.3±2.8mJy at FWHM 
widths of Aufwhm=366±124 and 529±75 kms^^ for the 
C0(J=2^1) and CO(J=3-^2) lines, respectively (see 
Tab.U). The line widths are consistent within the errors. 
A comparison of the line proflles may suggest more com- 
plex structure than a single Gaussian, but more sensi- 
tive observations are required to investigate this in more 
detail. The observed-frame peak velocities correspond 
to a median redshift of zco=2.0590±0.0003. We find 
velocity-integrated C0(J=2-J>1) and CO(J=3^2) line 
fluxes of /co=3.8±l.l and 9.1±1.1 Jy kms^i. The line 
brightness temperature ratio of r32 = 1.0±0.3 is consis- 
tent with optically thick, thermalized emission within 
the errors, consistent with w hat is found in other high-z 
quasars (e.g., Riechers et aL l2006( ). 

We derive a lensing-corrected CO line luminosity of 
L^o=1.3xl09(/iL/173)-iK kms-ipc^. Assuming a 
conversion factor of q;co=0-8 Mq (K kms^^ pc^)~^ from 
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Fig. 1.— CARMA CO(J=2->l) and CO(J=3^-2) maps {left) and spectra (right) toward B1938+666 (2=2.059). Contours arc plotted 
in steps of Scr (1(7=0.77 and 0.52 mjy beam~^ over 808 and 773 kms~^ (203 and 292 MHz)), starting at ±3(7 (no continuum emission 
subtracted). The beam sizes are indicated at the bottom of each map. The crosses indicate the (common) pointing center. The middle 
spectrum (histogram) shows the full covered frequency range of the redshift search from 74.1 to 115.3 GHz at a resolution of 62.496 MHz. 
The top and bottom spectra (histograms) of zoomed-in regions around the CO( J=2— ^1) and CO( J=3— >2) lines are shown at resolutions of 
20.832 and 10.416MHz (83 and 28 kms^^ for the CO J=2— ^-l and 3—^2 lines). The solid black curves are Gaussian fits to the spectra. 



L'qq to gas mass for ultra- l umino us infrared galaxies 
(ULIRGs; Downes & Solomon 'Tggs') , we find 
Mgas=1.0xlO^(;/XL/173)-i (aco/0.8)A/a_lsee Tab. O. 
In their analysis, Barvainis & Ivison (|2002D introduce 
a cutoff in /iL of a factor of 20, which may be more ap- 
propriate if the gas reservoir extends over a ^100 pc size 
region. Under this assumption, L'qq and Mgas would be 
~8x higher, but still at the low end of observed values 
at high z. 

3.1.3. New Results: Constraints on HCN, HCO+ , HNC, 
C2H, and CN Line Emission 

Based on the firm CO redshift, we went back to the 
spectrum to search for emission from dense molecular 
gas tracers. Fixing the redshift to z=2.0590 and all 
four lines to a common width (the width of the C2H 
line was corrected up by 2.6% to account for hyper- 
fine st ructure; see Fig. [21 and Riechers et al. I2007al . 
[2009a, for more details on the fitting procedure), we 
simultaneously fitted Gaussian line profiles at the red- 
shifted frequencies of HCN(J=3^2), HCO+ (7=3-5-2), 
HNC(J=3-J'2), and C2H(iV=3^2) emission0 adding 

4 The rest frequencies of the HCN, HCO+, and HNC(J=3^2) 
lines are i/rest=265.886180, 267.557619, and 271.981142 GHz. The 



one free parameter to fit the underlying continuum (see 
Fig. [3). We find S',.=3.1±l.l, 1.6±1.1, 1.9±1.1, and 
1.7±l.lmJy at a common At;FWHM=496±192 kms~^ 
for the HCN, HCO+, and HNC(J=3-l'2) and 
C2H(7V=3— >-2) lines, respectively. This corresponds to 
1.62±0.54, /hco+=0-86±0.45, /hnc=1-00±0.46, 
and /c2H=0.89±0.4(3 Jy kms~^, respectively. We also 
searched for CN emission, simultaneously fitting the 
CN(Af=3^2) and C0( J=3-^2) fines and the underlying 
continuum emission (see Fig. [S]). The width of the CN 
line was corrected up by 10% to account for hyperfine 
structure (see Fig. [1])(3 The fit yields 5'^=1.5±1.1 mJy. 
corresponding to /hcn=1-0±0.7 Jy kms^^. 

The HCN( J=3^2) line is marginally detected at 3.0a 
level, but we conservatively treat it as an upper limit 
in the following. The other four lines yield 1.4-2.2(7 
signals, and thus are not detected. Translating their 
fluxes to 3<T upper limtis yields /<1.4 (HCO+, HNC, 
C2H) and <2.1Jykms~i (CN). Given that the simul- 
taneously fitted HCN, HCO+, HNC, and C2H fines all 



C2H(7V=3— >2) line has six spectrally distinguishable hyperfine 
components at !/rost=262.004227-262. 208439 GHz (see Fig.^. 

^ The CN(Af=3— >-2) line has 19 spectrally distinguishable hyper- 
fine components at i/reat=339.446777-340.279166 GHz (see Fig.[2ll. 



CO in Lensed z>2 Quasars 



5 



85.4 



85.5 



Vobs [GHz] 
85.6 85.7 



85.8 



85.9 





A 




o _ 




1.6 - 


E 




> 

CO 


1.2 - 


> 




CO 


0.8 - 




0.4 




- 



Model spectrum of C2H(N=3-2) in B1 938+666 



261.2 




111.4 



2 Model spectrum of CN{N=3-2) in B1 938+666 



CO 

> 

CO 




339.2 339.6 



340 



340.4 340.8 341.2 



^rest 



[GHz] 



Fig. 2.— Model spectra of the C2H(Af=3^2) and CN(7V=3-i>2) 
emission toward B 1938+666. The relative intensities are co mputed 
in loc al thermodynamic equilibrium (see Riechers et al. |2007a| . 
I2009al , for details on the computation). The plus signs indicate 
all hyperfine components, for each of which the same line width as 
for the CO( J=3— >'2) line is assumed. The crosses and dashed lines 
indicate the peak intensities and profiles of the three main subcom- 
ponents (labeled "1" to "3"), summed over hyperfine components 
close in frequency. The vertical axes indicate the predicted inten- 
sity, normalized to the brightest subcomponent. The two brightest 
of these subcomponents for each line are indicated in Fig. |^ The 
solid lines are the profiles summed over all components. The thick 
dotted line indicates single Gaussian profiles fitted to the summed 
profiles, and have widths that are by ~2.6% (C2H) and ~10% (CN) 
larger than that of the individual components (and are used to fit 
the observations to minimize the number of fitting parameters). 

yield positive signals at 1.9-3.0(7 level, we attempted to 
stack their spectra. This yields a weighted average flux 
of 1.05±0.24 Jy kms"^, corresponding to a signal of 4.4(t 
significance. This stacked signal corresponds to 19%±5% 
of the C0( J=3— >2) luminosity (i.e., not corrected for line 
excitation). This is consistent with the median HCN/CO 
L' ratios (a measure of the dense gas fraction in galax- 
ies) i n nearby ULIR Gs and z>2 quasars (e.g., Gao et al. 
[2007t Riechers et aL l2007bD . 

3.1.4. New Results: Continuum Emission 

As shown in Fig. [TJ we detect continuum emis- 
sion toward B1938+666 over the entire observed wave- 
length range. Simultaneous fits to the line and 
continuum emission suggest continuum strengths of 
23.89±0.32 and 17.76±0.26mJy below the C0(J=2-^1) 
and CO(J=3^2) lines, and 23.15±0.15mJy below the 
HCN(J=3^2), HCO+(J=3^2), HNC(J=3^2), and 
C2H(A^=3-^2) lines. The strength and slope of the emis- 
sion is consistent with synchrotron emission from the 




Fig. 3.— A Search for HCN(J=3-i-2), HCO+(J=3-i>2), 
HNC(J=3^-2), C2H(Ar=3^-2), and CN(Af=3-s-2) emission toward 
B1938-I-666 (z=2.059). The spectra (histogram) are shown at the 
same resolution as in the middle panel of Fig. U Top: The solid 
black curve shows a simultaneous Gaussian fit to all four lines with 
a common line width (corrected for C2H hyperfine structure; see 
Fig. [2J and the underlying continuum emission, with the redshift 
fixed to 2=2.0590. The peak frequencies of the HCN(J=3^-2), 
HCO+(J=3^-2), and HNC(J=3->2) Unes as well as the main hy- 
perfine structure components of the C2H(Af=3— >2) line are indi- 
cated by the vertical lines. Bottom: The solid black curve shows a 
simultaneous Gaussian fit to both lines with a common line width 
(corrected for CN hyperfine structure; see Fig. [Jjl and the under- 
lying continuum emission. The peak frequencies of the main hy- 
perfine components of CN(Ar=3^-2) and that of the CO(J=3->-2) 
line are indicated by the vortical lines. 

radio-loud AGN. The spectrum may suggest a more com- 
plex continuum slope than a simple power law. How- 
ever, due to variability of the phase calibrator over the 
4 month course of the observations, there is a residual 
uncertainty in the flux calibration between the differ- 
ent frequency settings. We thus conclude that the con- 
tinuum slope is consistent with a power law within the 
errors. By averaging the line-free channels over the en- 
tire observed frequency range and fitting a 2-dimensional 
Gaussian to the u — v data, we find a 3.1 mm continuum 
flux of 20.84±0.14mJy. In Fig.|4l a map of the averaged 
continuum emission is shown. From the fit, we also find 
a continuum size of (1.48"±0.05") x (1.10"±0.05") at a 
position angle of 85° ±5°. Accounting for interferometric 
seeing (which is not part of the above errors), this is con- 
sistent with the 0.95" diameter of the Einstein ring, and 
the continuum emission to be dominated by the radio- 
loud AGN. 

3.1.5. Origin of the CO/Continuum Emission 

In Figure [5l an overlay of the C0( J=3— !>2) and un- 
derlying continuum emission on top of observed-frame 
1.6 /im continuum emission in B1938+666 is shown {Hub- 
ble Space Telescope (HST) image from King et al. 11991 
all HST images are data products from the Hubble 
Legacy Archive) The peak of the CO /continuum emis- 
sion is consistent with the brightest rest-frame opti- 
cal (523 nm) emission region along the Einstein ring. 
This suggests that the optical and (rest frame) (sub- 
)millimeter emission (observed frame 2.6-4.0 mm) are 
likely cospatial, with the optical emission being domi- 
nated by stellar light above the Balmer break. The lack 
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Fig. 4. — 3.1mm continuum emission in B13594-154 and B1938+666, averaged over all line-free data in the lower and upper sidebands. 
The beam sizes are indicated at the bottom. The crosses indicate the same position as in Figs. [T] and [7] Left: Contours are plotted in steps 
of 2(T (lo-=0.12 mjy beam~^ ), starting at ±3(t. Right: Contours are plotted in steps of SOcr (1(7=24 ^Jybeam"-'^), starting at ztSOcr. 
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Fig. 5. — Overlay of CO(J=3— >2) (contours) and 1.6 fim (rest-frame 523 nm) continuum emission {left; i75T/NICMOSl; image not 
cleaned; the background source is not detected at 555 and 814nm), and zoomed-in intensity niap of the 1.6 /^m emission {right; smoothed 
with a 0.1" Gaussian kernel) in B1938-f666. The cross indicates the same position as in Fig. [T] The source is lensed into a full Einstein 
ring, and the central spot is the foreground lensing galaxy. 



of pointlikc images along the Einstein ring suggest that 
the radio-loud AGN is optically obscured (or, alterna- 
tively, has low optical luminosity). The strength and 
spectral slope of the rest-frame submillimeter continuum 
emission make it unlikely that the AGN is much less 
magnified than its host galaxy. Thus, B1938-I-666 likely 
hosts a type-2 AGN along with its star burst. 

3.1.6. Absorption Line Searches 

Given the high 3 mm continuum flux of B1938-)-666, 
we have searched the spectrum for absorption lines from 
foreground sources. In Fig. [6l three spectral regions 
are shown that cover the absorption line frequencies of 



HCN(J=1^0), HCO+(J=1^0), and HNC(J=1^0) in 
the Milky Way (i.e., at z=0), and of HCN(J=2^1), 
HC0+(J=2^1), and HNC(J=2^1) at z=0.881, the 
redshift of the lensing galaxy. HCN, HCO+, or HNC 
are abundant molecules typically found in dense regions, 
and typically produce some of t he deepest a bsorption fea- 
tures (e.g., Wikhnd & Combes 119951 (l996[) . The spectra 
are normalized to an arbitrary flux scale, and re-binned 
to 10.416 MHz (35 and 33 kms""'^), to enable searches for 
narrow lines from compact molecular structures. At this 
resolution, the continuum emission is detected in every 
single velocity channel, but at relatively moderate signal- 
to-noise ratio. We do not detect any absorption lines. 
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TABLE 3 

Luminosities, gas masses and star-formation rates (SFRs). 



Source 


[Gpc] 




r b , c 
^FIR 

[1012 Lq] 


r/ b,d 
^CO 

[10l°K kms-ipc2] 


^fir/^CQ 
[Lq/K kms-lpc2] 


Mga.(H2)'= 

[IOIOMq] 


SFR* 
[Moyr-i] 


[Myr] 


B1938+666 


16.31 


173 


33/0.19 


(22.0±2.5)/0.13 


150 


0.10 


30 


35 


HE 0230-2130 


17.36 


14.5 


21/1.5 


(21.9±3.1)/1.5 


100 


1.2 


220 


55 


HE 1104-1805 


18.91 


10.8 


16/1.5 


(22.4±3.4)/2.1 


70 


1.7 


220 


75 


B1359+154 


28.36 


118 


11/0.093 


(7.7±1.2)/0.065 


140 


0.052 


14 


40 



^ Adopted from Barvainis & Ivison |20021 . 

^ Apparent luminosities (not corrected for lcnsing)/intrinsic luminosities (Icnsing— corrected). 
^ L-piR calculated assuming similar SED shapes as for the Cloverleaf quasar (WeiB et al. |200!j) . 

Extrapolated to CO( J— 1^0) luminosity by correcting for excitation, based on excitation modeling of high-js quasars (Riochcrs ot al. 12009 . l2009al) . 
Corrections applied: 0% for CO(J=2^1), 2% for CO(J=3^-2), 6% for CO(J=4^3). 



^CO(l-O) 



^ Assuming a conversion factor of a — 0.8 Mq /K kms ^pc^ from 
[1993) . 

^ Assuming (•Kcnnicutt ll998al : [T998H) : SFR[M0 yr~^]= 1.5 x 10"^" LfirILq], 
describes the dependence on the mass function of the stellar population, (5gB 
and not the AGN. 

^ Gas depletion timescales, defined as Tdcp — A^gas/SFR. 



to A/gas (H2) as appropriate for ULIRGs (see Downes & Solomon 

i.e., (5mfi5sb — 1-5 following the notation of Omont et al. 120011 : (5mf 
2;ives the fraction of Lpifj that is actually powered by the starburst 
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dv=35km/s 



[GHz] 



Lensing Galaxy: z=0.881 absorption lines 




HCN(J=l-a) 



HCO^(J=l-S) 



^v^33km/s 



94.5 

[GHz] 



Lansing Galaxy: z=0.881 absorption lines 




HNC(J.l-2) 



"«„ [GHz] 

Fig. 6. — Absorption line searches toward the continuum of 
B1938-I-666. The top panel shows a frequency region cov- 
ering the HCN(J=1^0), HCO+{J=1<-0), and HNC(J=l<-0) 
lines at z=0 (i.e., Galactic regions) at a velocity resolution of 
35 kms~l (10.416 MHz). The middle and bottom panels show fre- 
quency regions covering the HCN(J=2<— 1), HCO+( J=2<— 1), and 
HNC(J=2<-1) lines at 2=0.881, i.e., the redshift of the lensing 
galaxy, at a velocity resolution of 33 kms"-*^ (10.416 MHz). All 
spectra are normalized to an arbitrary flux scale. 

Assuming a velocity width that is at least comparable to 
the resolution of these spectra, this suggests that <20% 
of the continuum emission of B1938-I-666 are absorbed by 
Galactic or nearby extragalactic foreground structures, 
and <40% are absorbed by molecular gas (HCN, HCO+, 
or HNC) in the lensing galaxy. 

3.2. HE 0230-2130: The Southern Cloverleaf 
3.2.1. Previous Results 

HE 0230-2130 is an optically selected, quadruply 
lensed quasar at an optical redshift of Zopt=2.163±0.003 
with a maximum image separation of 2.15" (Wisotzki 



et al. 119991 Anguita et al. 120081), similar to the well- 
known Cloverleaf quasar (Magain et al. 1988^. It 
is lensed by two galaxies at almost identical redshift 
(zGi^ 0.523±0.001, ZG2=0.526±0.002; Eigenbrod et al. 
|2006[) . This 'Southern Cloverleaf is radio-quiet, but has 
a luminous dust bump, with 6*450 ;im=77±13mJy and 
'5'850pm=21.0±2.2mJy. These fluxes are ~35% of those 
found for the Cloverleaf, and consistent with the same 
spectral slope (Barvainis & Ivison 12002). Barvainis & 
Ivison find a model-based lensing magnification factor 
of /iL=14.5, which is similar to th at of the Cloverleaf 
(AiL=ll; Venturini & Solomon [2003). HE 0230-2130 also 
exhibits luminous emission from polycyclic aromatic hy- 
drocarbons (PAHs; Lutz et al. 2008). Together, these 
properties suggest the presence of a substantial molecu- 
lar gas reservoir. 

3.2.2. New Results 

We detect luminous, marginally resolved C0( J=3^2) 
emission toward HE 0230-2130 at >9(j significance 
(Fig. [T]) . From Gaussian fitting to the line profile, we find 
5*^=11. l±1.5mJy and Awfwhm=705±123 kms~^, cor- 
responding to /co=8.3±1.2 Jy kms~^(see Tab. [2]). We 
do not detect the underlying continuum down to a 2a 
limit of <1.5mJy. The CO line emission yields a host 
galaxy redshift of zco=2.1664±0.0005, which is consis- 
tent with the optical redshift of the quasar within the 
errors (dz=||zco~-Zopt||— 0.003). 

The CO line emission in HE 0230-2130 appears re- 
solved both spatially and in velocity. There is a clear 
difference between the CO line width at peak position 
(dashed spectrum in Fig. ^ and the spatially integrated 
line width. In Figure |8l a first moment map is shown 
(clipped below la in the individual velocity channels). 
Despite the moderate signal-to-noise ratio of this map, a 
clear velocity gradient is apparent. Accounting for beam 
convolution, the spatial and velocity structure is consis- 
tent with the lensed size of the source as defined by the 
maximum images separation of 2.15" within the uncer- 
tainties. 

Assuming a 2% correction for subthermal ex- 
citation of the C0(J=3— !>2) lincQ (based on 

^ This model-based, small correction factor is not statistically 
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Fig. 7.— CARMA CO(J=3^2) and CO(J=4^.3) maps (left) and spectra (right) toward HE0230-2130 (2=2.166), HE 1104-1805 
(z=2.322), and B1359+154 (x=3.240). Contours are plotted in steps of Itr (0.52, 0.65, 0.97, and 0.71 mjy beam-i over 1114, 450, 344, 
and 430 kms~^ (406.25, 156.25, 93.75, and 156.25 MHz)), starting at ±3(7, ±2cr, ±3(7, and ±3(7 (no continuum emission subtracted). The 
beam sizes are indicated at the bottom of each map. The crosses indicate the pointing centers. The spectra (histograms) are shown at a 
resolution of 31.25 MHz (86, 90, 115, and 86 kms"!). The solid black curves are Gaussian fits to the spectra. The dashed histogram in the 
top right panel indicates the peak position spectrum only. 
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RIGHT ASCENSION (J2000) 

Fig. 8. — First moment map of the CO(J=3— >2) emission in 
HE 0230-2130 shown in Fig. Velocity contours are in steps of 
At)=100 kms"^ 

CO line radiative transfer mod el ing o f liigli-2: 
quasars, e.g., Riechers et al. [200l l2009bD we de- 
rive L^o=1.5xlOiO(AiL/14.5)-iK kms-ipc2, and 
Mgas=1.2xlOiO(/iL/14.5)-i(aco/O.8)M0 (see Tah.^. 

3.3. HE 1104-1805 
3.3.1. Previous Results 

HE 1104-1805 is an optically selected double image 
quasar with a lens image separation of 3.19" at an 
optic al redshift of z t-=2 3192±0.0007 (Wisotzki et al. 
119951 : Sulentic et al. I2006D . It i s lens ed by a galaxy at 
ZG=0.729±0.001 (Lidman et al. |2000( ). It s host galaxy 
is lensed into an Einstein ring (Peng et al. 120061) . It has 
a supermassive black hole mass of Mbh=2.4x 10^ M©, 
and follows the relation between black hole mass and 
rest-frame i?-band luminosity of the host g alaxy (Mbh- 
Lr) for z>1.7 galaxies (Peng et al. I20060 Estimating 
its bolometric AGN lumin osity from its Mbh and B- 
band magnitude (Peng et al. l2006D . HE 1104-1805 follows 
the Lboi-ipiR relation for PC quasars, and is consis- 
tent with the values found for other high-z FIR-luminous 
quasars with in the uncertainties (e.g., Omont et al. 120031 : 
Wang et al. [20081: see Tab. El for Lfir). HE 1104-1805 
is radio-quiet, but has a luminous dust bump, with 
5*850 /jm=14.8±3.0mJy and 5i.3,nm=5.3±0.9mJy. It is 
magn ified by a factor of /zl=10.8 (Barvainis & Ivison 
I2002D. I t also exhibits luminous PAH emission (Lutz et 
al. l2008n . 

3.3.2. New Results 

We detect luminous, spatially resolved C0(J=3— 7-2) 
emission toward both lens images of HE 1104-1805 at 
>8ct and >4cr significance (Fig. [71). From Gaussian 
fitting to the line profile, we find S'i,=16.1±2.4mJy 
and Awfwhm=441±81 kms^^, corresponding to 
/co=7.5±1.2Jy kms^i (see Tab. ^. We do not de- 
tect the underlying continuum down to a 2a limit of 
<1.5mJy. The CO line emission yields a host galaxy 
redshift of zco=2.3221±0.0004, which is close to the 
optical redshift of the quasar (dz=||zco~-Zopt||— 0.003). 



Assuming a 2% correction for subthermal ex- 
citation of the C0(J=:3— >2) line, we derive 
L'co=2.1xl0i°(AiL/14.5)"iK kms^^pc^, and 
Mga,=1.7xl0iO(AiL/10.8)-i (aco/0.8)Mo (see Tab. [3]). 

3.3.3. Origin of the CO Emission 

In Figure [11 overlays of the C0(J=3— >2) emission on 
top of observed-frame 555 nm, 814nm, and 1.6 /im con- 
tinuum emission in HE 1104- 1805 are sho wn (H ST im- 
ages from Remy et al. Il998t Lehar et al. I2000D . The 
CO emission is clearly associated with the two (point- 
like) images of the AGN. The emission peaks are slightly 
spatially offset. This may indicate that the AGN does 
not reside in the center of the molecular gas reservoir. 
However, this small offset may as well be due to smaller- 
scale dynamics of the molecular gas that are averaged 
over in the CO emission line map. The CO images are 
consistent with being resolved along the Einstein radius, 
but observations at higher spatial resolution and sensi- 
tivity are desirable to study the structure and dynamics 
of the lensed gas reservoir in more detail. 

3.4. B1359+154 
3.4.1. Previous Results 

B1359-I-154 is a rare, sextuply lensed, radio-loud 
quasar at an optical redshift of 

•Zopt=3.235±0.002, with a m aximum image separation 
of 1.71" (Myers et al. |1999[) . It was selected in the 
radio, but in contrast to B1938-I-666, is an optically 
bright quasar. It is lensed by a compact group at zq— 1 
with three primary lensing galaxies of similar luminosi- 
ties, which are situated on the vertices of a triangle, 
separated by ^^0.7" (photometric/lensing redshifts are 
ZGi=1.35±0 16, ZG2=0.88±0.06, and zg3=0.94±0.07; 
Rusin et al. |2001[ ). Very compact (3-9 mas lensed sizes 
along their major axes), fiat spectrum (af y~-0.3 from 
1.7 to 5 GHz) radio cores are associated with all six im- 
ages, and the brightest three images show radio jets ex- 
tending out to tens of mas (lensed) scales (Rusin et al. 
1200 1|) . B1359-I-154 has a luminous dust bump, with 
5*450 /im=39±10mJy and 5850Mm=ll-5±1.9mJy. It is 
magn ified by a factor of /xl=118 (Barvainis & Ivison 
[2002). 

3.4.2. New Results: Line Emission 

We detect luminous CO(J=3^2) and CO(J=4^3) 
emission toward B1359+154 (Fig. [7]). From Gaus- 
sian fitting to the line profiles, we find 5,^=5. 6±1. 7 
and 10.0±1.6mJy, and Aufwhm=198±92 and 
237±47kms-i for' the CO(J=3^2) and CO(J=4^3) 
lines, respectively. The line widths are consistent 
within the errors. This corresponds to /co=l-2±0.4 
and 2.5±0.4Jy kms-i (see Tab. [2]) . The CO line 
emission yields a median host galaxy redshift of 
zco=3.2399±0.0003, which is close to the optical 
redshift of the quasar (dz=||zco--Zopt||— 0.005). The 
line brightness temperature ratio of r43=1.2±0.5 is 
consistent with optically thick, thermalized emission 
within the errors. 

Assuming a 6% correction for subthermal excitation of 
the (higher signal-to-noise ratio) C0(J=4— >-3) line, we 
derive L'co=6.5xl08 (/LtL/118)-i K kms'^pc^, and 
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Fig. 9. — Overlays of CO(J=3— >2) (contours) and 555 nm, 814 nm, and 1.6 /xm (rest-frame 167, 245, and 482 nm) continuum emission 
{H3T/WFPC2 and NICMOS2, left to right; images not cleaned) in HE 1104-1805. The cross indicates the same position as in Fig. [7] The 
two brightest spots are the lens images of the background quasar, and the fainter spot in the middle is the foreground lensing galaxy. 



o 15°13'25' 




35:5 

RIGHT ASCENSION (J2000) 

Fig. 10. — Overlays of CO(J=4— >3) (contours) and 555 nm, 814 nm, and 1.6 ^m (rest-frame 131, 192, and 377 nm) continuum emission 
(HST/WFPC2 and NICMOS2, left to right; images not cleaned) in B1359-I-154. The cross indicates the same position as in Fig. [7] The six 
spots visible in all panels are the lens images of the background quasar, and the three additional spots in the central region that become 
visible at longer wavelengths are the foreground lensing galaxies. 

for interferometric seeing (which is not part of the above 
errors), this is consistent with the 1.7" maximum image 
separation, and the continuum emission to be dominated 
by the radio-loud AGN. 



Mgas=5.2xlG« 



(AiL/118)-i(aco/O.8)M0 (see Tab. El), 
comparable to what is found in the lensed Ly-break 
galaxy MS 1512-cB58 (Riechers et al. I2010b( ). Under 
the assumption of /zl=20, L^q -^gas would be <6x 
higher, but still at the low end of observed values at high 



3.4.3. New Results: Continuum Emission 

As shown in Fig. [71 we detect continuum emission to- 
ward B1359+154 under the CO lines. Simultaneous fits 
to the line and continuum emission suggest continuum 
strengths of 3.97±0.35 and 2.11±0.32mJy below the 
CO(J=3^2) and CO(J=4^-3) lines. The strength and 
slope of the emission is consistent with synchrotron emis- 
sion from the radio-loud AGN. By averaging the line- free 
channels over the observed frequency range (LSB-I-USB 
for both lines) and fitting a 2-dimensional Gaussian to 
the u — V data, we find a 3.1mm continuum flux of 
3.61±G.18mJy. In Fig. HI a map of the averaged con- 
tinuum emission is shown. From the fit, we also find 
a continuum size of (3.14"±2.48") x (0.88"±0.72") at a 
position angle of 289°±19°, suggesting that the contin- 
uum emission is marginally resolved at best. Accounting 



3.4.4. Origin of the CO/Continuum Emission 

In FigurefTUl overlays of the C0( J=4— )-3) and underly- 
ing continuum emission on top of observed-frame 555 nm, 
814 nm, and 1.6 /im continuum emission in B135 9-I-154 
are shown {HST images from Munoz et al. Il998| ). The 
peak of the CO/continuum emission is consistent with 
the centroid of the six lensed images of the source. Ob- 
servations at (at least) 5x higher resolution are desirable 
to study the distribution and dynamics of the molecular 
gas in this complex system. 

4. DISCUSSION 

4.1. The 'Star Formation Law' for AGN Host Galaxies 

Due to high lensing magnification factors, our targets 
are among the intrinsically faintest high- 2: sources ob- 
served to date. Their Lpm/ico ratios lie at the low end 
of the distribution of high-z quasars, but are consistent 
with the observed range. 
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Fig. 11. — Comparison of CO line luminosity with FIR luminosity for samples of low- 2 : spiral and starburst galaxies and (U)LIRGs (Gao 
& Solomon 2004), low-z HE and PG QSOs (Evans et al. 2001, 2006; Scovillc ct a l. |2003 : Bertram et al.[2007), and high-z S MGs (Solomon 
& Vanden Bout 2005; Tacconi et al. 2006; lono et al. 2006; Frayer ct al. 200i|, [2Q3; Schinnerer et al. 2008; Daddi et al. 2009a, '20093; 
Knudsen et al. 2009; WeiB et al. 2009; Bothwell et al. 2010; Carilli ct al. 2010; Riech ers et al. l2010al : Harris ct al. 2010) and QSOs (Solomon 
& Va nden Bout i2005i ; R iecher s et al.^2006, 2009a; Carilli et al. ,2007t Maiolino et al. l2007l : Willott et al. i2007i : Aravcna et al. .2008; Coppin 
et al. 120081 ; Wang et al. 120101 ). The new lensed high-z QSOs in this paper are shown as double circles. For comparison, the green circles 
assume /iL=20 for B1359+154 and B1938-f666. All data are corrected for gravitational lensing. The solid line is a fit to all quasar samples 
(black symbols), corresponding to logLpij)^ = (1.33±0.05) log LJ-,q-(1.09±0.50). The dashed line is a fit to all sources without a dominant 
AGN (magenta symbols), corresponding to logLpia = (1.47±0.06) log L'(2q-(2.73±0.54). A fit to all data yields LpiR = (1.40±0.04) 
logL^O"(l-92±0.42). 



To investigate the properties of our targets in more 
detail, we here compare their L'qq and £fir to other 
CO-detected quasars at low and high redshift, and low 
and high redshift galaxies without luminous AGN see 
Fig. [TT]) . The L'^Q-L-pm, relation is a variant of the 
spatia l ly integ rated 'star formation law' (e.g., Kennicutt 
Il998at Il998b( ) that only relies on observed quantities, as 
it is independent of quantities like the conversion factor 
aco and the assumed stellar initial mass function (see, 
e.g., Daddi et al. 2010a; Genzel et al. 2010, for recent in- 
vestigations of the Mgas-SFR relation). At high redshift, 
we restrict the analysis to a comparison of quasars (sum- 
marized in Table HJ and submillimeter galaxies (SMGs), 
as these samples cover a comparable range in LfirQ 
Many SMGs are known to host AGN, but these AGN are 

* The high-2 massive, gas-rich star-forming galaxies identified 



substantially less luminous than those found in quasars. 
At low z, the galaxies without luminous AGN are repre- 
sented by a revised version of the Gao & Solomon (200|) 
sample, including luminous and ultra-luminous infrared 
galaxies ((U)LIRGs). This currently is the best-studied, 
best calibrated low-z galaxy sample for the spatially in- 
tegrated galaxy properties studied here. We corrected all 
high-J CO line-based luminosities for excitation, based 
on models for high-z quasars with well-sampled CO line 
ladders (Riechers et al. 2006; 2009b), corresponding to 
corrections of 0%, 2%, 6%, 12%, and 25% for the CO 
J=2^1, 3->2, 4^3, 5^4, and 6->5 lines. A fit to 
aU data yields ipiR = (1.40±0.04) logL'co-(l-92±0.42), 



by Daddi et al. l l200g. I2010bl ) and Tacconi et al. l l20l(D have com- 
parable L'qq, but substantially lower L^ib./ L'qq, than these more 
extreme systems. 
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Fig. 12. — Lpiji/ L'^Q, a measure of star formation efficiency, as a function of Lpia (left) and redshift {middle), and as histograms for 
the low-2 and high-2 sample s (r ight). The ratio is shown on a linear scale to emphasize the scatter of the ^qq^^fiR relation. The same 
samples are shown as in Fig. 1111 



yieldi ng the c anonical power law slop e of 1 .4 (e.g., Ken- 
nicutt Il998at see also Riechers et al. |2006[) . Fitting all 
quasar samples yields logipiR = (1.33±0.05) logL^Q- 
(1.09±0.50), i.e., a slightly lower power law slope. A fit 
to the low-z galaxy sample and the SMGs yields logLpiR 
= (1.47±0.06) logi^o-(2.73±0.54), i.e., a slightly higher 
power law slope. 

The high-z quasars and SMG samples show compa- 
rable Lfir/^qq ratios and span a comparable range 
(see also Fig. [I2l where the ratio of both quantitites is 
shown on a linear scale to emphasize the scatter of in- 
dividual sources). On average, SMGs have somewhat 
higher ratios (by ~28% at Lfir=10"'^'^ L©), which is 
likely due to the fact that the typically C0(J=3— i>2)- 
based L^q are under-corrected with the 2% correction 
factor for quasars. Recent studies of C0(J=1— J'O) and 
C0( J=2— >1) emission in SMGs suggest that a correction 
factor of 30%-40% may be more appropria te for a sub- 
stantial fraction of SMGs (e. g., Iviso n et al. 120101 : Carilli 
et al. 120101: Riechers et al. I2010al: Harris et al. 120101 : 
Frayer et al. I2010D . which would make the LpiR/i'co 
ratios for quasars and SMGs even more similar on aver- 
age. However, studies of the CO excitation in SMGs also 
show that the emission seen in low- J lines may be par- 
tially associated with a low-excitation gas component, 
i.e., material that likely is not concentrated in the re- 
gio ns wh ere the massive s tarbur sts take place (Carilli et 
al. 120101: Riechers et al. l2010aD . This is not the case 
in quasars, which can typically be modeled with a sin- 
gle, highly excited gas co nipone nt down to the low - J CO 
lines (e.g., Riechers et al. 120061 : Weifi et al l2007f ). Stud- 
ies of C0( J=l— ?>0) emission in larger samples of high-z 
quasars and SMGs are required to investigate these dif- 
ferences, which directly affect the physical interpretation 
of the Lqq-Lfir relation. In fact, observations of CO 
line ladders down to CO(J=1^'0) may be a powerful 
tool to distinguish different high-z galaxy populations 
ba sed on their gas prop erties alone (e.g., Da nnerbau er et 
al. 120091 : Aravena et al. I2010t Riechers et al. I2010b[ ). 

Another important aspect is that we here study Lfir 
(i.e., rest-frame 42. 5-122. 5/im emission), not Lir (8- 
1000 /im). This has a relatively small impact for galax- 
ies without luminous AGN, but it makes substantial 



difference for quasars. Quasars typically show both a 
warm dust bump with temperatures of rdust=30-60K 
(i.e., comparable to ULIRGs and SMGs), and a hot dust 
bump with Tdust >100 K (likely associated with the AGN, 
e.g., Weifi et al. [20031 Beelen et ah [20061 ). This yields a 
Lin/Lpm ratio of typically a factor of a few. Thus, a 
L'qq-Lib, relation for quasars would likely be biased to- 
ward the AGN properties. On the other hand, our find- 
ings for the Lqq-Lfir relation appear to suggest that, 
indeed, the FIR luminosity in high-z quasars is domi- 
nated by dust-reprocessed emission from young stars in 
the host galaxy, rather than the AGN. 

One caveat to the present interpretation of the L'qq- 
LpiR relation at high z is the possibility of differen- 
tial gravitational lensing effects between the CO and 
FIR emission for the lensed galaxies in both samples. 
Given the overall agreement between lensed and unlensed 
sources within the samples and the fact that they follow 
the Lqo~^fir relation, such effects appear to be minor, 
but only few lensed systems have been studied at suf- 
ficient spatial resolution in molecular gas and FIR con- 
tinum emission to quanitfy this in more detail. 

The difference in Lfir/^co ratios between low-z 
quasar samples and low-z galaxies without luminous 
AGN appear to be larger than those between high-z 
quasars and SMGs. With the current sample sizes, these 
differences are only moderately statistically significant. 
Given the less massive dust and gas reservoirs in the 
lower luminosity systems, such a trend may however not 
be surprising. In such systems, the AGN may be sub- 
stantially more efficient in heating larger fractions of the 
warm dust than in the typically very gas-rich high-z sys- 
tems, yielding elevated LFiR/i'co ratios. High spatial 
resolution observations of the CO and FIR continuum 
in nearby quasar host galaxies are key to investigate 
this issue in more detail. An alternative interpretation 
would be that the AGN fraction rises with Lfir toward 
ULIRGs and SMGs, which could decrease the difference 
between FIR-luminous quasars and other FIR-luminous 
galaxy populations. We consider this explanation less 
compelling, as there is no clear relation between AGN 
luminosity and FIR luminosity for galaxies of compara- 
ble masses. 
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TABLE 4 

Molecular Line and FIR continuum properties of high-z quasars. 



Source 


2CO 


ML 


CO line 


[Jy kms~-^] 


logio(L^P,)'''' 
[K kms pc^] 


logio(LFIR)'' 

[Lq] 


ref. 


comments 


Q0957+561 


1.4141 


1.6 


2-s>l 


1.2+0.1 


10.30 


12.94 


1 


double lens 


3C318 


1.571 




2^-1 


1.19+0.22 


10.59 


13.00 


2 




COSBOll 


1.8275 




2^1 


1.33 


10.76 


12.93 


3 




B1938+666 


2.0590 


173 


3^2 


9.1+1.1 


9.10 


11.28 


4 


Einstein ring 


HS 1002+4400 


2.1015 




3^2 


1.7+0.3 


10.63 


13.08 


5 




HE 0230-2130 


2.1664 


14.5 


3—5-2 


8.3+1.2 


10.18 


12.16 


4 


quadruple lens 


RX J1249-0559 


2.2470 




3—5-2 


1.3+0.4 


10.56 


12.85 


5 




F10214+4724 


2.2858 


17 


3—5-2 


3.40+0.19 


9.76 


12.55 


6 


arc lens 


HE 1104-1805 


2.3221 


10.8 


3^2 


7.5+1.2 


10.32 


12.17 


4 


double lens 


J1543+5359 


2.3698 




3^2 


1.0+0.2 


10.49 


13.00 


5 




HS 1611+4719 


2.3961 




3^2 


1.7+0.3 


10.73 


13.08 


5 




J13122+12381 


2.5564 




3^2 


0.4+0.1 


10.15 


12.48 


5 




Cloverleaf 


2.55784 


11 


3^2 


13.2+0.2 


10.63 


12.73 


7 


quadruple lens 


J1409+5628 


2.5832 




3—5-2 


2.3+0.2 


10.92 


13.28 


8 




MG 0414+0534 


2.639 


27 


3^-2 


2.6 


9.55 


12.07 


9 


quadruple lens 


RX J0911+0551 


2.796 


22 


3^2 


2.9+1.1 


9.73 


12.37 


10 


quadruple lens 


J04135+10277 


2.846 


1.3 


3—5-2 


5.4+1.3 


11.25 


13.38 


10 


cluster lens 


MG 0751+2716 


3.1999 


16 


3->2 


4.6+0.5 


10.17 


12.49 


11 


quadruple lens 


B1359+154 


3.2399 


118 


4^3 


2.5+0.4 


8.81 


10.97 


4 


sextuple lens 


APM 08279+5255 


3.9118 


4.2 


1^0 


0.168+0.015 


10.40 


13.68 


12 


triple lens 


PSS J2322+1944 


4.1192 


5.3 


1^0 


0.155+0.013 


10.30 


12.86 


13 


Einstein ring 


BRI 1335-0417 


4.4074 




2->l 


0.43+0.02 


10.91 


13.45 


14 




BRI 0952-0115 


4.4337 


4 


5^4 


0.91+0.11 


9.89 


12.38 


15 


double lens 


BR 1202-0725 


4.6949 




1^-0 


0.120+0.010 


11.00 


13.78 


13 


two sources 


J0338+0021 


5.0267 




5^4 


0.73+0.09 


10.48 


13.18 


16 




J0927+2001 


5.7722 




5^-4 


0.44+0.07 


10.35 


12.93 


17 




J1044-0125 


5.7824 




6^5 


0.21+0.04 


9.92 


12.72 


18 




J0840+5624 


5.8437 




5^4 


0.60+0.07 


10.50 


12.86 


18 




J1425+3254 


5.8918 




6^5 


0.59+0.11 


10.38 


12.73 


18 




J1335+3533 


5.9012 




6^5 


0.53+0.07 


10.34 


12.74 


18 




J2054-0005 


6.0379 




6->5 


0.34+0.07 


10.16 


12.74 


18 




J1048+4637 


6.2274 




6->5 


0.39+0.08 


10.24 


12.77 


18 




J1148+5251 


6.4189 




3^2 


0.20+0.02 


10.48 


13.13 


19 





[6] 



this work; [5] Coppinet 
[1 1] All oin et al. I2007t 
al. 120071 : [17] Carilli et aj. 120071 : [18] Wang et aL I20I 



Ao et al. 1200a 
Riechers et aL I2009d 



Referen ces. — - [1] Krips et al. 120051: [2] Willott ot al. 2007;_[3] Aravcna ct al. .2008,; [4] 
Woifi ct al. I200a [8l Beelen et al. 120041 : [9] B arvainis et al. 119981 : [1 0] Ha inline et al. I200l 
Ricchors ct al. 1200^ [1 4] Riechers et al. 120081 : [15] Guilloteau et al. 119991 : [16] Maiolino et 
[19] Walter ct al. book 

Note. — Estimates are typically based on the lowcst-J CO line observed, if applicable. 

^ Corrected for gravitational Icnsing. 

^ Corrected for line excitation where applicable. Corrections applied: 0% for CO(J=2^1), 2% for CO(J=3-5-2), 6% for CO(J=4-5-3), 12% for 
CO(J=5^4), and 25% for CO(,7=6-5-5). 



4.2. Star Formation Rates and Gas Depletion 
Timescales 

If the LpiR in our targets are indeed dominated by 
star formation, we can derive their star-formation rates 
(SFRs). Assuming SFR[M ^ vr-i]= 1.5xl0-i° ipraii©] 
(e.g., Kennicutt Il998al: ll998bD . we find SFRs 
of 3O(AiL/173)-iM0yr-i, 220 M^yi-^ 
220(/ZL/10.8)-iMQyr-\ and 14 (^l/IIS)-^ yr"! 
for B1938+666, HE0230-2130, HE 1104-1805, and 
B1359+154, respectively (Tab.©. 

The minimum times for which the starbursts can be 
maintained at their current rates are given by the gas 
depletion timescales. We find Tdep=Mgas/SFR=35, 55, 
75, and 40Myr for B1938+666, HE 0230-2130, HE 1104- 
1805, and B1359-(-154, respectively. This is comparable 
to what is found in SMGs and other high-z QSOs (e.g., 
Greve et al. 120051 : Riechers et al. 2008; 2010a,). 

4.3. Systemic Redshifts 

We find that the CO redshifts of the three sources 
with spectroscopic redshifts from optical observations are 



close to 



'opt • 



However, in all cases, the CO emission is 



ically CivA1549) optical emission lines from the AGN. 
The median redshift difference of the three sources with 
Zopt is dz=0.0045±0.0015, or 360±116 kms-^EI For a 
large sample of Sl oan D igital Sky Survey (SDSS) quasars, 
Richards et al. ()2002[ ) have found that low-ionization 
lines such as Mgii(A2798) are commonly redshifted rel- 
ative to high-ionization broad lines such as Civ(A1549) 
(median of 824±511 kms~^), which they interpret to be 
due to orientation effects related to the AGN accretion 
disk (or disk wind) . Even though the orientation of the 
accretion disk and host galaxy need not be aligned, refer- 
encing the redshifts of broad AGN lines to the systemic 
redshifts provided by CO observations of the host galax- 
ies may help to better disentangle orientation effects and 
dynamics of high redshift galaxies. 

5. SUMMARY 

We have detected C0(J=2-^1), CO(J=3^2), and 
C0( J=4— J-3) emission in four lensed quasar host galaxies 
at redshifts of z>2 with CARMA (a total of six lines). 

^ For consistency, we here adopt 2Lya=2.3172 for HE 1104- 
1805 as measured by Smette et al. 1)19951 ) instead of the 



redshifted relative to broad Lya and high ionization (typ- H/3+[Oiii] (A5007) redshift used above. 
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To date, these are the highest redshift CO detections re- 
ported with CARMA. From our analysis, we obtain the 
fohowing key resuhs: 

Facihtating the large, >7.4GHz bandwidth 
(LSB+USB) of the new CARMA correlator, we 
have executed the first successful 'blind' CO redshift 
search with an interferometer, demonstrating the fea- 
sibility of such investigations with only few frequency 
setups. Based on the CO redshift, we also used the 
broad frequency range covered by these observations to 
search for emission from dense molecular gas tracers 
(HON, HCO+, HNC, C2H, and CN), which yielded only 
a statistical detection after stacking - as expected at the 
depth of this search. We also searched for HCN, HCO"*", 
and HNC absorption toward foreground sources, but no 
lines are detected along the line of sight. 

We spatially resolve the CO emission in the two radio- 
quiet quasars in the sample. In the wide separation dou- 
ble lens HE 1104-1805 (z=2.322), we individually detect 
the emission toward both images. In all cases the CO 
emission is consistent with the brightest emission regions 
detected at rest-frame optical wavelengths. 

We detect strong continuum emission toward the two 
radio-loud sources in the sample. The spectral slope of 
the emission is consistent with synchrotron emission as- 
sociated with the AGN in both cases. 

We derive lensing- and excitation-corrected CO line 
luminosities of 0.65-21 x 10^ K kms~^ pc^. From the CO 
luminosities, we determine molecular gas masses of 0.52- 
ITxIO^Mq. These values are at the lower end of those 
observed in other high-z quasars, showing the advantage 
of observing lensed distant galaxies to probe down to 
lower gas masses. 

Combining our targets with literature samples, we find 
no significant difference in the L^yrJ L'qq ratios between 
high-z quasars and SMGs. We find tentative evidence 
that nearby quasars with low LpiR show an excess in 



LfyrI L'qq relative to systems without a luminous AGN 
at similar redshift and with comparable L'qq. This 
may indicate that AGN heating of the warm dust is 
more efficient than in the more gas-rich and massive 
FIR-luminous quasars detected at high z, where the 
dust heating in the FIR appears to be largely dominated 
by young stars. 

The observations presented herein lay the foundation 
for future studies with the Atacama Large (sub-) Mil- 
limeter Array (ALMA), which will routinely observe un- 
lensed galaxies at comparable and fainter intrinsic CO 
luminosities. The 'blind' CO redshift search technique 
presented here will be particularly valuable for so-called 
'Deep Field' investigations with ALMA, which will allow, 
for the first time, to construct an unbiased CO luminos- 
ity function out to high z by directly selecting galaxies 
through their CO content. This is an important step to- 
ward constraining the gas mass history of the universe, 
a critical piece in our understanding of galaxy evolution 
throughout cosmic times. 
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